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Abstract. Dark matter annihilation in Galactic substructure produces dif- 
fuse gamma-ray emission of remarkably constant intensity across the sky, and in 
general this signal dominates over the smooth halo signal at angles greater than 
a few tens of degrees from the Galactic Center. The large-scale isotropy of the 
emission from substructure suggests that it may be difficult to extract this Galac- 
tic dark matter signal from the extragalactic gamma-ray background. I show 
that dark matter substructure induces characteristic small-scale anisotropies in 
the diffuse emission which may provide a robust means of distinguishing this 
component. I discuss predictions for the angular power spectrum of the diffuse 
emission from various extragalactic source classes as well as from Galactic dark 
matter, and show that the energy dependence of the angular power spectrum of 
the total measured emission could be used to confidently identify a gamma-ray 
signal from Galactic dark matter. 



1. Introduction 



The observation of gamma-ray emission from the annihilation of dark matter 
particles could provide a means of indirectly detecting dark matter, constraining 
its properties, and mapping its spatial distribution. T he Large Area Teles cope 
aboard the Fermi Gamma-ray Space Telescope {Fermi, lAtwood et al.l [20091 ) will 
measure for the first time the all-sky gamma-ray emission at energies from a few 
tens of GeV to a few hundred GeV, an energy band which is especially promising 
for indirect dark matter searches. 

A generic prediction of cold dark matter models is the persistence of sub- 
structure within dark matter halos. Depending on the details of the particle 
physics model, the smallest dark matter halos for rn in the early universe with 
masses of lO'^^ _ jg-e fe.g.. iGreen et al.ll2005l : IProfumo et "alll2006l ^. and 
many of these dense halos are expected to survive to the present day. This 
prediction is consistent with the results of numerical simulations, which find an 
abundance of subhalos within a Galaxy-sized dark matter halo. Simulations also 
find that the radial distribution of subhalos is less centrally concentrated than 
the smooth dark matter of the host halo (e. g., Gao et al.P 2004: Diemand et al 



120071 : iKuhlen et aDl2007l : iMadau et all l2008l : ISoringel et al..,2008al ). and conse- 
quently the angular distribution of Galacti c subhalos ori the s ky as observed 
frorn our position appears fairly isotropic ( Kuhlen et al.l [20081 : ISiegal-GaskinsI 
[20081 : [Springel et al.[[2008b[ ). as shown in the left panel of Figure [H 
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The rate of dark matter annihilation scales with the square of the density, 
so nearby regions of high density, such as the Galactic Center and subhalos, are 
natural targets for indirect searches. While the Galactic Center is almost cer- 
tainly the brightest point source of gamma-rays from dark matter, substantial 
astrophysical foregrounds make extracting a dark matter signal from the Galac- 
tic Center a formidable challenge. Subhalos, on the other hand, are largely free of 
such foregrounds, and although most subhalos in the halo of t he Milky Way are 
not likely to be reso lved individually as point sources by Fermi ( Fieri et al.ll200^ : 
Kuhlen et al.ll2008l ). collectively substructure is expected to generate a large flux 
of diffuse gamma-rays, exc eeding the smooth halo flux for angles > 10 degrees 



from the Galactic Center (iPieri et al.1 l2008l : iKuhlen et all l2008l : ISoringel et al 
2008bl : iFornasa et al.l |2009| ) . The right panel of Figure [1] shows a realization 



of the signal from substructure as observed by an experiment with the angular 
resolution of Fermi above ~ 10 GeV. 

The 'isotropic' diffuse gamma-ray background thus may include contribu- 
tions from Galactic dark matter substructure and extragalactic dark matter, as 
well as from unresolved members of known extragalactic source classes such as 
blazars. Since most of the measured dark matter flux may be hidden in dif- 
fuse emission, it is worthwhile to explore methods for robustly extracting this 
signal. Recognizing that different source classes have different anisotropy prop- 
erties, several recent studies have considered anisotropics in diffuse emission on 
small angular scales as a way of statistically identifying contributions from dis- 
tinct source classes. Angular powe r spectra of the gamma,-ray emission frorn 
several extragalact i c source cl asses fAndo &i Komatsu' '2006'; 'Ando et al.' '2007|; 



Miniati et al.ll2007l : ICuoco et a l. 2008; Taoso et al. 2 009; F o rnasa_et_al. 20 09) as 



well as from Galactic dark matter (jSiegal-GaskinsI 1200^ : iFornasa et al.1 l2009l : 
Anddl2009l ') have been predicted. 



In this paper, I review recent theoretical work related to anisotropy anal- 
ysis of gamma - ray d ata, however the content presen t ed is based primarily on 
Siegal-Gaski"^ (|2008l 1 and lSieeal-Gaskins k Pavhdoul (|2009 '). In Section 2 I dis- 



cuss predictions for the angular power spectra of individual source classes and 
their relevance to the angular power spectrum of the total measured emission. 
I consider the anisotropy energy spectrum, or the energy dependence of the 
total angular power spectrum, in Section 3 and show that the combination of 
anisotropy and spectral information can be a powerful tool for disentangling con- 
tributions to the diffuse emission from multiple or changing source populations. 
I summarize my conclusions in Section 4. 



2. Angular power spectra of gamma-ray source populations 

The angular power spectrum provides a statistical measure of fluctuations as 
a function of angular scale. Here I consider the angular power spectrum Ci of 
intensity fluctuations 51{ij)) = {I{ip) — (/))/(/), where is the intensity in 
the direction ij). Note that the fluctuations are expressed in units of the mean 
map intensity, hence 51 is dimensionless. The angular power spectrum is given 
by = ( la^mP), where are determined by expanding 51{'ijj) in spherical 
harmonics, 5I{'4)) = Y,i^m^f--mY£m{ip)- 
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Figure 1. Gamma-ray intensity from dark matter annihilation in Galactic 
substructure. The right panel shows the same map as the left panel, but 
smoothed with a Gaussian beam of width 0.1°, cor responding to Fermi 's 
angular resolution above ^ 10 GeV. Figures are from lSiegal-GaskinsI (|200l . 



For an unclustered distribution of point sources (i.e., a distribution for which 
the angular position of a given source is uncorrelated with the position of the 
other sources), the angular power spectrum is noise-Hke. In this case, Ci takes 
the same value at all i, and is inversely proportional to the number density of 
sources per solid angle. Blazars are point sources and only weakly clustered, 
hence their angul ar power spectrum is well-approximated as noise-like (see, e.g., 
Ando et al. 20071 ). Most Galactic subhalos are point-like due to the extremely 



small angular size of their emission region, and from our position do not ap- 
pear strongly clustered, so Galactic dark matter substructure also produces a 
fairly noise-like power spectrum. However, subhalos which are sufficiently mas- 
sive or nearby are extended sources, and their contribution results in relatively 
less angular power at high multipoles (corresponding to small angular scales) 
compared to a pure noise-like source distribution. In general, the amplitude of 
the angular power spectrum from Galactic substructure is predicted to be much 
greater than that from blazars. This can be understood by noting that as a 
cosmological source class, blazars have a much higher number density per solid 
angle than Galactic substructure, and thus produce less angular power. 

Figure [2] shows the angular power spectrum from Galactic substructure 
for subhalo populations which extend to a minimum subhalo mass M^[^ = 10 
Mq and Mmin = 10^ Mq, and for the smooth dark matter halo. In the small 
Mrain scenario, the angular power spectrum is quite noise-like, reflecting that 
most of the emission is from subhalos that are effectively point sources. For 
the large M^[^ case the angular power spectrum is flatter because more of the 
emission originates from extended sources. The overall normalization of the 
power spectrum is higher in the Mmin = 10^ Mq case since the mean map 
intensity is smaller, and thus the same intensity fluctuation is a larger fractional 
deviation from the mean than in the M^[^ = 10 Mq case. Due to its small 
contribution to the total intensity, the smooth halo contributes minimally to the 
total angular power. It is important to keep in mind that because the angular 
power spectra shown are calculated in units of the mean intensity of each source 
class, to compare the contributions of these source classes to the total angular 
power spectrum the individual power spectra must be scaled by the fractional 
contribution of the intensity of that source class squared. 
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Figure 2. The angular power spectrum from Galactic dark matter substruc- 
ture for Afniin = 10 M(7), Mniin = 10 ^ , and the smooth dark matter halo. 
Figure is from iSieeal-Gaski^ (l2008h . 

In general, the intensity and angular power spectrum of the emission from 
Galactic dark matter substructure exceeds that from extragalactic dark mat- 
ter, and among known extragalactic source classes, unresolved blazars are the 
dominant contributor to both the intensity and angular power spectrum at en- 
ergies > a few hundred MeV. Here I consider only these two source classes, as 
they are expected to be the most relevant for analysis of the diffuse gamma-ray 
background with Fermi. 

3. The anisotropy energy spectrum 

In this section I present two example scenarios that demonstrate that the energy 
dependence of the angular power spectrum of the total measured emission could 
be used to confidently infer the presence of multiple contributing source pop- 
ulations. In these examples I assume that the total measured isotropic diffuse 
emission is comprised exclusively of emission from unresolved blazars and from 
Galactic dark matter substructure. For cosmological source populations such as 
blazars, at energies above a few tens of GeV interactions with the UV, optical, 
and IR backgrounds (the extragalactic background light, EBL) lead to an expo- 
nential cutoff in the observed intensity energy spectrum, so for the blazar spectra 
EBL attenuation is approximated here by applying an energy-dependent cut-off, 
and assuming for simplicity that all blazars reside at a single redshift. The dark 
matter intensity e nergy spectrum is calculated using the approximation given in 
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The angular power spectrum of the total diffuse emission from extragalactic 
sources Cf^ and Galactic dark matter substructure Cf^ is given by 

where /eg and /dm are the (energy-dependent) fractions of the total emission 
from extragalactic sources and Galactic dark matter substructure, respectively. 
Extragalactic sources and Galactic dark matter substructure are uncorrelated, 
so here the cross-correlation term (^^GxDM _ q_ evident from Eq. [H the 
amplitude of the total angular power spectrum at a given energy depends on 
the relative contributions of each source class to the measured emission at that 
energy. For clarity, I use the term intensity energy spectrum to describe the 
energy dependence of the photon intensity. 

The anisotropy energy spectrum is defined as the value of the angular power 
spectrum at a fixed multipole i as a function of energy. The angular power 
spectrum of a single source class in which all members have the same intensity 
energy spectrum is independent of energy because the angular distribution of 
the sources on the sky is identical at all energies. However, contributions to 
the diffuse emission from multiple populations with different intensity energy 
spectra result in a modulation of the anisotropy energy spectrum, as the relative 
contribution of each population changes with energy. 

Figures [3] and H] present the intensity energy spectra and corresponding 
anisotropy energy spectra at £ = 100 for two example scenarios. In each scenario 
the intensity energy spectrum of the total measured diffuse emission is consistent 
with that of an unresolved blazar population plus a dark matter component, 
but also with that of an unresolved blazar population with slightly different 
spectral properties without a dark matter component. In these cases spectral 
information alone is insufficient to identify the dark matter signal, however, the 
anisotropy energy spectrum can be used to robustly identify the presence of a 
dark matter contribution. In both examples the observed 'total' intensity energy 
spectra and anisotropy energy spectra are the sum of the dark matter and the 
EBL-attenuated reference blazar contributions. The intensity energy spectrum 
of an alternative blazar model is also shown which, with EBL attenuation, could 
plausibly produce the observed total intensity energy spectrum without a dark 
matter component. 

The scenario in Figure [3] assumes a dark matter mass of 700 GeV, and de- 
spite dark matter dominating the total intensity above ~ 20 GeV, its contribu- 
tion is not clearly identifiable in the intensity energy spectrum. In the anisotropy 
energy spectrum, on the other hand, the dark matter contribution produces in a 
significant feature. In Figured] the dark matter mass is assumed to be 80 GeV, 
and here dark matter never dominates the intensity energy spectrum, but still 
produces a strong feature in the anisotropy energy spectrum. In both examples 
the anisotropy energy spectrum enables the blazar plus dark matter scenario to 
be distinguished from a blazar-only scenario: if unresolved blazars were the sole 
source of the isotropic diffuse emission, the anisotropy energy spectrum would 
be constant in energy, but the presence of a dark matter contribution results 
in a modulation in the anisotropy energy spectrum, unambiguously indicating a 
change in contributing source populations with energy. 
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In these figures the error bars indicate the 1 — a statistical uncertainty in 
the measured angular power spectrum for an all-sky observation time of 5 years 
with Fermi, assuming a fraction /g^y = 0.75 of the sky over which the analysis 
is performed. The 1 — a uncertainty is given by 



here corresponds to Fermi's angular resolution. The data is assumed to be 
averaged over multipole bins of width A£ = 100. The noise power spectrum 
Cn is the sum of the Poisson noise of the signal and the background, Cn = 
(4vr/sky/A's)(l + (A^b/-^s)), where A'g and iVb are the number of signal and 
background photons observed. Unresolved blazars and Galactic dark matter 
substructure are considered signal, while any other observed emission, such as 
Galactic sources other than substructure, are noise. The error bars shown gen- 
erously assume N]^/Ns = 10, i.e., an order of magnitude more background than 
signal photons. 

For the anisotropy energy spectrum at i = 100, the error bars rapidly 
increase at low energies (E < 1 GeV) due to the degradation of Fermi's angular 
resolution below this energy, and at high energies due to low photon statistics. 
At intermediate energies the noise term in Eq. [2](Cn/VF/) is quite small, and in 
this regime a modulation in the measured anisotropy energy spectrum can be 
easily identified. 

It is notable that in the optimistic case that the baseline extragalactic 
anisotropy level is observable, it may be possible to extract the shape of the 
dark matter annihilation spectrum via Eq. [H even in cases where the dark mat- 
ter signal cannot be disentangled from other contributions in the intensity energy 
spectrum. 

4. Summary 

Anisotropies in diffuse gamma-ray emission can be a powerful diagnostic of the 
properties of contributing source populations. In particular, indirect searches for 
dark matter may be able to use fluctuations on small angular scales to statisti- 
cally detect dark matter substructure in the halo of the Galaxy even if individual 
subhalos are not detectable. 

The robustness of anisotropy analysis can be improved by examining the 
energy dependence of anisotropies, since the angular power spectrum of the 
total measured emission varies with energy according to the energy-dependent 
contribution of each source population. The examples presented in Section 3 
demonstrate that a contribution from Galactic dark matter substructure to the 
measured diffuse emission leads to a departure of the anisotropy energy spectrum 
from an energy-invariant quantity. For plausible scenarios, this modulation is 
shown to be detectable in Fermi data with high confidence. 

Acknowledgments. It is a pleasure to thank my collaborator, Vasiliki 
Pavlidou, for her invaluable contributions to the work presented. 
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Figure 3. Example intensity energy spectrum and corresponding anisotropy 
energy spectrum assuming a dar k matt er particle mass of 700 GeV. Figure is 
from .Siegal-Gaskins fc Pavlidoui (|2009f ). 
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Figure 4. S ame as Figure Bl for a dark matt er particle mass of 80 GeV. 
Figure is from iSieeal-Gaskins fc Pavlidoui (|2009l) . 



8 



J. M. Siegal-Gaskins 



References 

Ando, S. 2009. larXiv:0903.4685 l 

Ando, S., & Komatsu, E. 2006, Phys. Rev. D, 73, 023521 

Ando, S., Komatsu, E., Narumoto, T., & Totani, T. 2007, Phys. Rev. D, 75, 063519 
Atwood, W. B., at al. 2009, ApJ, 697, 1071 

Bergstrom, L., Edsjo, J., & Ullio, P. 2001, Physical Review Letters, 87, 251301 
Cuoco, A., Brandbyge, J., Hannestad, S., Haugb0lle, T., & Miele, G. 2008, 

Phys. Rev. D, 77, 123518 

Diemand, J., Kuhlen, M., & Madau, P. 2007, ApJ, 667, 859 

Fornasa, M., Pieri, L., Bertone, G., & Branchini, E. 2009, arXiv:090T.2921 

Gao, L., White, S. D. M., Jenkins, A., Stoehr, F., & Springel, V. 2004, MNRAS, 355, 

819 

Green, A. M., Hofmann, S., & Schwarz, D. J. 2005, Journal of Cosmology and Astro- 
Particle Physics, 8, 3 
Kuhlen, M., Diemand, J., & Madau, P. 2007, ApJ, 671, 1135 
Kuhlen, M., Diemand, J., & Madau, P. 2008, ApJ, 686, 262 
Madau, P., Diemand, J., & Kuhlen, M. 2008, ApJ, 679, 1260 
Miniati, F., Koushiappas, S. M., & Di Matteo, T. 2007, ApJ, 667, LI 
Pieri, L., Bertone, G., & Branchini, E. 2008, MNRAS, 384, 1627 

Profumo, S., Sigurdson, K., & Kamionkowski, M. 2006, Physical Review Letters, 97, 
031301 

Siegal-Gaskins, J. M. 2008, Journal of Cosmology and Astro-Particle Physics, 10, 40 
Siegal-Gaskins, J. M., & Pavlidou, V. 2009. ■ar Xiv:0901.3776l 
Springel, V., et al. 2008a, MNRAS, 391, 1685 
Springel, V., et al. 2008b, Nature, 456, 73 

Taoso, M., Ando, S., Bertone, G., & Profumo, S. 2009, Phys. Rev. D, 79, 043521 



